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TECHNICAL MEMORANDUM X- 64865 


EFFECT OF GRAVITATIONAL AND AERODYNAMIC 
TORQUES ON A RIGID SKYLAB-TYPE SATELLITE 

I. INTRODUCTION 


The behavior of a Stylab-type satellite under the influence of gravita- 
tional and aerodynamic torques is investigated. Only the simple case of free 
molecular flow of uniform velocity and completely diffuse reflection (i.e., the 
atmospheric particle transfers its total relative momentum to the satellite at 
impact) is considered. The aerodynamic torque on the satellite is expanded 
into a Fourier series with respect to the two angles characterizing the direction 
of the atmospheric velocity vector relative to the satellite. Making use of some 
invariance properties and symmetries of the satellite, and neglecting some 
relatively small effects of the torque, these Fourier series can be simplified; 

i.e., it can be shown that many terms vanish. The resulting expressions, 
suitably truncated, coincide with those given by Nurre [1] , which he derived 
by numerical methods. Applying the standard methods of linear stability to the 
linearized equations of motion, it is found that the satellite does not have a 
stable equilibrium position. The presentation comprises sufficient details to 
make it a useful aid in similar studies. 


II, ASSUMPTIONS 


Throughout this study, the following assumptions are made: 

1. The satellite is a rigid body. 

2. The satellite orbits the central body in a circular orbit; the central 
body is represented by a point mass, orbital and rotational motion of the satel- 
lite are not coupled, and only the principal term of the gravitational torque is 
taken into account. 

3. The atmosphere is resting with respect to a frame with inertial 
directions (i.e., the primary and its atmosphere do not rotate). 



4. Only the case of completely diffuse reflection (accommodation 
coefficient equals 1) is studied, i.e., the atmospheric particle is at rest with 
respect to the satellite immediately after hitting the satellite; in other words, 
the atmospheric particle transfers its entire relative linear momentum to the 
satellite. 


III. EQUATIONS OF MOTION 


Three coordinate frames will be used; an inertial frame, an orbital 
frame, and a satellite-fixed body frame. All coordinate systems are right- 
handed cartesian systems. 


A. Inertial Frame Y 

The inertial frame Y has its origin at the center of mass of the satel- 
lite, and the directions of the coordinate axes are space fixed. The Y^-axis 
and Y^-axis are in the plane of motion, and the Y 2 -axis is in the direction of 
the orbital angular momentum vector. The inertial frame will not be used 
explicitly, but serves as a basic reference and starting point to develop 
coordinate transformations and equations of motion. 


B. Orbital Frame F 

The orbital frame i has its origin at the center of mass of the satellite, 
and it rotates uniformly with respect to the inertial frame with angular speed 
S2. The Ij-iixis coincides with the line from the central body to the satellite, 
the Ij-axis is in the direction of the satellite's velocity, and the coin- 

cides with the Y 2 -axis. We have 


1 t 0 

-sS2 t^ 

— — 


0 1 

1 

0 

Y, = 

a 

y t 0 

cQ t j 




where is the constant angular velocity of the orbital frame with respect to 
the inertial frame, represented in the orbital frame (Fig. l). 
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Figure 1. Inertial and orbital coordinate systems. 

C. Body Frame 

The body frame is fixed in the satellite; we will use the following four 
coordinate systems fixed in the body. 

1. ^-Coordinate System . This is a general body fixed system, with 
its origin at the satellite’s center of mass. Occasionally we will specialize 
it to utilize certain properties of the satellite, such as symmetries, etc. 

Many of our equations wUl be written in this system. The rotation from the 
orbital |"-system to the T}-system is represented as a product of three plane 
rotations: rotate through with fixed 1-axis, then through with fixed new 
2-axis, then through ip^ with fixed new 3-axis. Thus 
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I ^2 ^3 ^3 ^3 ^3 “ ^3 ^ 

^2 ®3 ^1 ^3 ^ ^1 ^2 ®3 ^3 ^1 ®2 ^3 


\ 


Si Cj 


OJ = 


! C2 Cg !/li + Sg lp2 
-Cg Sg i^i + Cg i^2 
\ S2 ipi 


Cl Cj 


I 


\ 


h j 


( 2 ) 


( 3 ) 


where we have written Cg for c^2» co is the angular velocity of the tj- 
system with respect to the orbital |‘-system, represented in the ^-system. 

2 . X-Coordinate System , The origin of the X-coordinate system is at 
the satellite's center of mass; the directions of the axes are chosen such that 
J^e inertia matr^ is diagonal (principal axes system). The rotation from the 
I -system to the X-system and its matrix ® are defined as for the t} system, 
with 0^1 , » ^3 replacing ^2^ h' 

3 « Z-Coordinate System . The origin of the Z-coordinate system is 
at the satellite’s center of mass; the directions of the axes are chosen so as to 
£xhibit geometrical symmetries, etc., of the satellite. The rotation from the 
^—system to the Z— system and its matrix $ are defined as for the tj— system, 
with 01, 02 > replacing 0i, 0g. 

^-Coordinate System . The axes of the ^-coordinate system are 
parallel to those of the 'Z-system; its origin is chosen so that the coordinates 
adapt to the geometry of the satellite. 

* 

^ If Z is th6 vector from the origin of the Z- system to the origin of the 
^-system, then 

— ^ 

S = z - Z . (4^ 
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Figure 2 illustrates the various coordinate systems and their relations. 



Figure 2. Diagram of coordinate systems and their relations. 

The equation of rotational motion of the satellite, referred to a body 
frame, reads 


lof + cu*x Iw*'- i = 0 


(5) 


where 

I = (constant) moment of inertia matrix, 

ok 

o) = angular velocity of the body frame with respect to the inertial 
frame (and represented in the body frame), 

i = i + I = gravitational torque plus aerodynamic torque, 
g a 

For the subsequent discussions, we shall represent all quantities in the body 
frame; we then have 

■/ 



CO = cu + , ( 6 ) 

where, as stated before. 
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a; - angular velocity of the body frame with respect to the orbital 
frame, represented in the body frame, 

— * 

^ = tpn is the angular velocity of the orbital frame with respect 

to the inertial frame, represented in the orbital frame; for rep- 
resentation in the body frame it is multiplied by the transforma- 
tion matrix ip. 

It follows that 


— ^ ^ 

CJ = Ct) 4- ^12 + 


JL* 

and with = 0 in our case. 


u? =" CO - CO X lj}il 


or 


— * ^ 

CO = CO - CO X fi 


We thus arrive at the equation of motion 


I (uj - cu X S2) + (cu + n) X I(u) = 0 

S a 


( 8 ) 


where 


^ = ipQ, = 


^22 

yp32 j 


= o 


/Cl Sg + S, Sg CgX 
Cl C3 - Si S2 S3 
\ - St Cg / 


O) is given in equation (3), and 


( 9 ) 
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IV. GRAVITATIONAL TORQUE 

For a satellite orbiting the central body in a circular orbit with constant 
angular speed the principal term, of the gravitational torque is 

Ig = 7 X l7 , ( 11) 

where r is the unit vector from the ce_ntral body to the satellite. For the 
orbital frame, denoting the vector by r*, we obviously have 



hence 
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The atmosphere is considered as a continuous medium which transfers 
its entire relative momentum to the surface of the satellite at impact. 

A, Differential Force 

We calculate the force exerted by the atmosphere on a surface element 
of the satellite in the body frame. The basic equation from which we start is 

dpT = vdm , (16) 
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where 


ju = (linear) momentum, 


V = velocity of atmosphere with respect to the satellite, 


m = mass of atmosphere. 

The following fairly obvious sequence of equations ( Fig. 3), 


dm = p dV , 

(17) 

dV = da dh , 

(18) 

dh = ds COSO! , 

(19) 

cosa = e e , 

(20) 

n V 

ds = V dt , 

(21) 


where 


p = density of atmosphere, 


V = volume of atmosphere. 


da = surface element of satellite, 
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~ normal unit vector of surface element, pointing toward the 
"outside," 

= v/v = unit vector of v. 


yields upon substitution an equation for d^. 



Figure 3. Atmospheric force on surface element. 


_ Obviously, only those parts of the surface will contribute to the momen- 
tum n that satisfy the following two conditions: 

1, The surface part is not shaded against the atmospheric stream by 
another surface part. 

2. COSO > 0, i.e., 0. ( 22 ) 


Combining these results, we find for the momentum transferred to the surface 
element da 


d^i = pv^ (_e e ) e da dt 
n V V 



e < 0 

V 


(23) 
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and for the force df = djj/dt 


df = pv^(-e e)eda,ee<0. (24) 

n V V n V 


For later use, we note that cosq: da = -e e da is the projection of the 

n V 

surface element da onto a plane orthogonal to the velocity vector v. 


B. Differential Torque 


From the force df we get the differential torque di with respect to the 

P 

point P as 


df = (R-P)xdf 
P 


or 

di = pv^ ( -e e ) ( R - P) X e da , e e o , ( 25 ) 

p nv 'v nv '' 

where R is the location vector of the surface element da. 

The total force f and torque i are calculated by integrating equations 

(24) and (25) over the permissible surface area; i.e., the area that satisfies 
the conditions (22). 


VI. FOURIER EXPANSION OF AERODYNAMIC FORCE AND TORQUE 


Except for very simply shaped bodies, the aerodynamic force and torque 
as given by (the integrated form of) equations (24) and (25) are complicated 
functions of the body shape and the direction of the velocity vector v of the 
atmosphere (via the integration over the body surface). It results, however. 
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that both force and torque are, for a fixed given body, periodic functions of the 
two angles describing the direction of the atmospheric velocity vector. 


We shall subsequently establish the Fourier series representing these 
periodic functions and in the next part specialize it to a Skylab-type satellite. 
Consider the atmospheric force and torque in the body-fixed coordinate sys- 
tem 77; then, T) = Ip ^ . In the orbital frame, the atmospheric velocity vector 


V 


I 


-V 



( 26 ) 


and for the 7 )-coordinate system 



Comparing this with the representation 


( 27 ) 


lea c/?\ 
cq; s/3 

\sa J 
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of a unit vector in spherical coordinates where o! is the latitude and /3 is the 
longitude, we that we can interpret ^2 latitude and ^3 as negative 
longitude in the 77-coordinate system. 

In general, the area of the satellite exposed to impact of atmospheric 
particles is a discontinuous fuhction of the direction of the atmospheric veloc- 
ity vector; i.e., of the components C2C3, -C2 S3, of e^, and thus of the 
angles ^2 ^3 (think, for example, of a planar part of the satellite that 

changes from not being exposed to being exposed to impact instantly at certain 
values of ^2 ^-^d Nevertheless, the projection of the area onto a plane 

orthogonal to the atmospheric velocity vector and the projection's boundary 
are continuous functions ( of the body shape and) of the components of e^. 

As mentioned previously, this projection of the surface element da is 


da* = -e e da (e e < O) 
n v ' n v 


( 28 ) 


hence 


force f (^2> ^3) = pv^ e J da* ( 29 ) 

projected 

area 


and 


torque £ (i^2» h) ~ f (R-^xe da* (30) 

^ projected ^ 

area 

= pv^ / Jr da*\ X e -pv^ P x e Jda* , ( 31 ) 
(projected! projected 

\ area / area 


where the seco^ term on the right of equation ( 3 l) is just - P x f . All inte- 
grals and thus f and are continuous functions of tp2 and ip^; to repeat, the 
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angles enter the functions £ and i only through the components of e in the 

P V 

combinations c ^2 ^^ 2 - We can now write the following Fourier 

series; 


00 

= s 

^,1^ = 0 


(A, 


fXV 


C/X ^2 ^^^3 


+ D 
ixv 


Spt^2 


(32) 


i = 



00 

Z (E^J^ cm!/'2 c;^^3 + 


SjU^2 


+ H 

fJtl^ 


Slllp2 SVlpi) 


(33) 


or, in easily understood abbreviation, 


f = Sfcc + 2 f cs + 2 fsc+ Sfss 


i = 2^ cc + S i cs + 2i sc + S ^ ss 


( 34 ) 


Equations ( 34 ) state that each component of f and J is of the indicated form, 
where, for example, Efcs stands for 2 (appropriate coefficient for f-term)* 
cp ^2 s^'^ 3 , and all sums are taken over the appropriate values of the indices, 
as later stated in each case. 

By various invariance considerations, some general and some special- 
ized to a Slq^lab-type satellite, we will show that”most’’ of the coefficients in 
equations ( 32) and ( 33) actually vanish. 

By previous remarks the Fourier series must remain unchanged under 
all transformations^f the angles ^2 ^3 ihat leave the components 0^2 c^ 3 > 

c ^2 sil >2 6 ^ unchanged. The set of equations 
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Ci^'2* = Cip2 Ci^3 

Clp2* Slp^* = Clp2 Si/ig 


( 35 ) 


Slp2* = S02 
has the two solutions 


(Trivial) I: 1 P 2 * = ^2 . h* = ^3 

II: ^ 3 * = ^ + ^3 


(36) 


hence we must have 


f (^2. h) = f(’T - ipz> ^ *h) 

^(i/'2, h) = '^(’T - ^2» 7T + Ips) 


(37) 


By a routine process of con^aring coefficients we find the reduced series 

H+v even fj,+ v even v odd ti+v odd 

I = - E isc + Z'^®® 

M+ V even ju+ u even /x+ v odd |u+ u odd 

The second consideration applies rigorously to the force F only. It 

is obvious that the projected area fda+ is the same for e and -e ; therefore 

v _ v 

[see equation (30)] the substitution implies 7 — -f The equations 


(38) 

(39) 
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( 40 ) 


Ctp2* o4>2* = -Ci/-2 Clp^ 
Cip 2 * = -C1P2 S1P2 


have the two solutions 


:: il)^* = -i^2 » 

i/lj* = 3T + i^3 

: 1P2* = IT + ^2 

. h* = h 


( 41 ) 


both these solutions lead to the reduction of equation ( 38 ) . to 

^ + Zf«^ Z^^^^ • (42) 

M odd pL odd /x odd /x odd 

t/ odd vodd. IX even ixeven 

VII. SPECIALIZATION TO A SKYLAB-TYPE SATELLITE 


We assume now that the satellite is symmetric with respect to the rjj, 
173-plane. This implies that for a surface element da at tjj, ^3 with 

~ ^®nl ^n2 ®n3^ there is a surface element at ij^, -t) 2, 173 with e^'^ = 

~^n2 ^n3^ same area. Write f as a sum of integrals over 

the "right” (rj2 < 0 ) and the ”lefr (r?2^ O) parts R and L of the satellite: 


IS 




Under the reflection 

(^1 V2 nz) - (r?i -tJ2 Vz) , e^2 ^3) ^ Kl ~%2 ®v 3 ^ 

corresponding to the tjj, 173-plane symmetry, the integrals change according 
to 

M 

f — /|eG,-e„e_ + e^e„jj-e„ da 
^ nl vl o2 v2 n3 v3 I v2 1 

WJ 

and ( 45 ) 

l%i\ 

\3%s' -^v2r‘* V 

L R / 

Ws/ 

Hence, under this reflection the force f is changed according to 

= (fifzfs) -(fi-fzfa) . (46) 
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as we would expect intuitively. Similarly, we find that the torque I is 
changed according to 


i - (i j ^2 -^3) (--^1 -^2 -^3) • (47) 


(observe that for equation (47 ) to hold we assume that the reference point P 
lies in the symmetry plane 772 =0.) 

We now use the just discussed symmetry for a further reduction of the 
Fourier series ( 39) and ( 42). The transformations of the angles il>, and 
corresponding to the transformation (44) are solutions of 


C^2* = C1P2 C1P2 

CIP2* Slp2* = -CJp2 Si^3 

si^2* = s^2 ; 

the two solutions are 


I; ^2* = h . i/'s* = -i/'s 

II: ^2* = Tr-ip 2 > i/'s* = 


(49) 


Taking into account the change of f and i according to (46) and (47), both 
solutions (49) lead to the same further reduced Fourier series 




fs = ^ fee + 
M odd 
u odd 


Z fsc 

M odd 
V even 


(50) 
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( 51 ) 


^2 = JfCS + Jfss 
H odd fi odd 
u odd u even 

■^1 ’-^ 3 “ ■*■ Z '^®® 

fjL+ u even fi+ v odd 


(52) 


^2 = Z^^^^ ^ Z 

M+veven /x+j/odd 


(53) 


To achieve another reduction of the force terms, we can argue as 
follows. Under the transformation 4>2~- -^2 (reflection on the r}^, ? 72 -plane) 
the projected areas remain unchanged (if we neglect shading effects, which 
are considered to be rjelatively small), since every part of the satellite is 
(assumed to be) symmetric with respect to some plane parallel to the tJj, 
/? 2 -plane (Fig. 4). Hence, if we change 


- T 

e 


V 


^^vl ^v2 ^v3^ 


(e e -e ) 
^ vl v2 . v3^ 


(54) 


which corresponds to 


f = ( fi f 2 fs) — ( fi f 2 -fa) ( 55) 


(see, for example, the foregoing discussion of the 'n^, ijj-plane symmetiy), 
and it follows that fj and f 2 are (approximately) even functions of ^ 2 » and f 3 
is (approximately) an odd function of Applied to equations (50) and ( 51), 
this leads to 


fi = Z 

fx odd 
V odd 


(56) 
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Figure 4. 77 -coordinate system with respect to satellite. 


f2 = Z fcs 
M odd 
u odd 


(57) 


fs = Z 

ju odd 
V even 


(58) 


To further reduce the Fourier series for the torque i, we have to use 
the specific geometiy of the satellite and must restrict ourselves to only 
approximately valid arguments and conclusions. Thus, we assume that the 
satellite is composed of a main cylindrical part 1, of circular cross section 
and with orthogonal planar end panels 2a and 2b, and of rectangular flat panels 
3a, 3b, and 5. The part 4, connecting the cylindrical mantle 1 and panel 5, 
will be neglected (Figs. 4 and 5). 
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Figure 5. Shape and parts, of idealized satellite. 

For the following considerations, the origin of the coordinate system 
is at the center of the main cylinder 1, and the 773 -axis is its centerline. The 
panel 2 = 2 a + 2b is parallel to the r/i, 772 -plane, and the panels 3 = 3a + 3b 
and 5 are parallel to the 172 , 773 -plane. Both panels 3 and 5 are symmetric 
with respect to the t 7 j, 773 -plane (assumed basic symmetry). [We argue here 
in an elementary geometrical and mechanical manner. A more precise deri- 
vation can easily be supplied by the method used to derive equation (46).] 
Shading effects are neglected. 

First, we proceed as before in a general f^hion^we get, however, 
results only for ig. Consider the transformation -e^. It is easily seen 

that the torque £ with respect to the origin of the panels 3 = 3a + 3b and 5 

transforms like T— -7. For the mantle 1 and panel 2 = 2 a + 2 b, the torque 
transforms like 7—7 (observe, for example, that for e^^ > 0 only 2 a, but 

not 2b, which is shaded by 1 and 2a, experiences torque). Because of sym- 
metry, both for 1 and 2 , we have ^3 = 0 ; thus 

ig s J -^cs+ J £ss . (59) 

H odd lA odd 

V odd V even 

Under the transformation ^2 ^ i*e., (e^^, e^^ e^^) — (e^^ e ^2 “% 3 )> 

both panels 3 and 5 we see that fg — ij (also — ij, which we cannot use); 
thus, fg must be an even function of ^ 2 * we get 

ig s 2 ics . (60) 

uodd 
p odd 
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A direct calculation of the torques on the various parts and addition, 
with shading effects neglected as before, yields considerably sharper results. 
Applying formula (31) for the torque 7, it can be shown that the torque of the 
complete cylinder, mantle 1, and panels 2a and 2b is always zero.’ Thus, we 
need to consider only the panels 3 = 3a + 3b and 5. Both 3 and 5 are sym- 
metric with respect to the plane ~ therefore, the vector q from the 
origin to the geometric center of 3 or 5 is 


q 



(61) 


The torqu^pf each panel with respect to its geometric center vanishes, and 
the force f can easily be computed from equation ( 29) to be 


f = pv^ A I e I e 

vl V 


(62) 


where A is the total area of the panel. Thus, the torque i on each panel is 
of the form 


-^3 % 2 \ 

"vl-‘il%3 • (®3) 

^v2 / 

We see that the total torque with respect to the geometric center of the 
cylindrical main body 1 + 2a + 2b on the satellite is of the same form. 

Generally, the Fourier series for Icoso- 1 Is a cosine series containing 
only even multiples of o; ; using this fact, we find 


i = pv2 A I e I q X e = pv2 A 1 e I 
vl v vl 


-^3 ~ ! C ^2 1 1 1 C^2 = 

E (coeff) cp ^2 

E(coeff) svi ^3 


p odd 

p odd 


= E^cs 
p odd 
V odd 


(64) 
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(65) 


! C 1 p 2 !l Clp.^ I (q 3 Clp 2 Clp 2 - Qi Slp^) = 


^icc + 2-^sc 
^ odd ju odd 

V odd V even 


We assume further that the center of mass of the satellite lies in the plane ^ 
172 = 0 geometric coordinate system (the specialized 17 - 

system used last) and the X-system (principal axes system) are parallel. If 
p is the vector from the origin of the X-system to the origin of the 77 -system, 
we get for the total torque i with respect to the center of mass 

C *0*111* 


r, rr, " P ^ ^ = -<!2 + PS ^1 " Pi ^3 

G * V/ *1X1* 



( 66 ) 


and using f|, f 2 ,_and fj from expressions (56) through (58), we see that the 

components of i have the same structure [expressions (64) and (65)] 

c.o.m. ^ 


as those of i. 


The expressions ( 56) through (58) for the force components and (64) 
and ( 65) for the torque components have, suitably truncated, the form given 
by Nurre [ l|/p. 1047, equations 5 1 (there, in the last term, read c 5®2 for s3(S^). 


Now we use the angles and for the form of the final expressions for 

7 and T, we take that given by Nurre [ 1 ] [for the notation, see expressions 
(32) and (33)]: 




f 2 




Am ^^2 ^^3 ■^131 ^^2 AgjI C 3 i ^2 ^^3 

B 112 Ci?2 ® ®132 ® *^2 Bgj2 c3i?2 ®332 C3i?2 S3t?j 

C103 ®<^2 ^123 + C303 s 3 i ?2 ^323 sSi?2 C2i?3 


>(67) 
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•^1 = Fill S.i^3 + Fi 3 i Ct ?2 s 3J^3 + Fgu C3,?2 6 i?g + F 331 C3i?2 S3i?3 | 

^2 ~ F112 C-1^2 Ci?3 + Ej 32 Ci? 2 cSt^g + Egjj c 3 i ?2 Ct?3 + E332 C 3 i ?2 

■*■ G 102 Si >2 ■•■ G 122 S ??2 C 2 i ?3 + G 3 Q 2 S3i?2 ■*■ G 322 S3i?2 c 2 i >3 

■s 

■^3 Fji 3 C^2 Si?3 + Fi33 Ci?2 SSiJg + F313 CSt^g SiJg + F333 C3i?2 SStJg j 


. ( 68 ) 


VIII. EQUILIBRIUM 


The equilibria of a system x = g(x) of differential equations are defined 
to be_^e special solutions with all derivatives set to zero; i.e., the solutions 
of g(^ — 0. The equilibria of the equation of motion (8) are the solutions of 


f2xIS2-^-f = 0 
g a 


where is given by ( ll) and ( 15). by (68) , and Q by (9). 


With 


I = 



we have 


22 ^32 


(i 3 -y« 

nxm . n? I 

JI 2 “ ^ 1 ) ^12 ^22 


(69) 


(70) 


(71) 
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The vector equation ( 69 ) represents a system of three scalar equations for 
the three quantities 1?^, general there exist finitely many 

solutions mod 2 tt. 

In our case, the solution of most interest is i?i = 1^3 = 0 ; the first and 
third components of ( 69 ) vanish identically, and the second is of the form (A, 
B, . . . , E are constant coefficients) 


ACi?2 Bs,?2 Cc3i>2 + Ds3i? 2 + Ec,>2 ® 


( 72 ) 


which can be solved for i?2; denote a solution of interest by i?* • 

IX. LINEARIZED EQUATIONS OF MOTION 


For the investigation of the linear stability of the previously described 
equilibrium = 1J3 = 0, " ‘^2 » we set 


i?i - 6 ^ , i?2 - 1?^ + 62 * ^3 “ ®3 


( 73 ) 


and linearize the equations of motion (S) with respect to the 6. and their 
time derivatives. Using ^ 


c(i?2 + ^2) = Ci^l - 62 s,?2 = c* - 62 s* 
s(^f + 62) - + 62 = s+ + 630* 


( 74 ) 


we easily derive 


s Ii + I 3 
6 i + J_^ 


I2 s* 


p. 


Ii + Iq - I9 1 


63 - 4 ^ 

c* 


h - h 


61 -(Ml + s^Llii s*) 63 = 0 


( 75 ) 
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(c*^ s*2) _ Mj 62 = 0 

^2 


63 + s*6i + — 2 — 1 c*6i + 4 ^ s* 5 i + ^_il(l + 3 s*2 ) _ Mgl 60 

^3 h L I 3 J 


where Mj, M2, and are certain constants. 


X. LINEAR STABILITY 


Equation ( 76) for 63 is not coupled to the other two equations; its 
stability behavior is determined by the coefficient of 5 ,. The solution is 
unstable if 


^2 


for > 0 in ( 78) the solution can be represented by trigonometric functions. 

Of more interest is the system ( 75), (77); it is of the form 

A6 + B6 + C6 = 0 , (79) 

where A, B, and C are 2 by 2 matrices and 6 ^ = 82). Introduce 3=6; then 

6 - 3 = 0 


A3 + C6 + B3 =0 
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or in matrix form, 



( 81 ) 


where U is the 2 by 2 unit matrix. Since 


u oyi /U 0 \ 

.0 a/ ’ \0 A'V 


we get the standard form 





(82) 


The coefficient of the cubic term in the characteristic polynomiai of ( 82) is 
the trace of the coefficient matrix, which is here the trace of A~^B; since 


A'‘ = 



B 


Ii + I3 - I2 








(83) 


we conclude that this trace vanishes. Moreover, this coefficient is the sum 
of the four roots of the characteristic polynomial, and since it vanishes, not 
all roots can have a negative real part (the sum of the real parts must vanish) 
and the solutions cannot be stable. 
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XI. CONCLUSIONS 


The discussions and results of Sections IX and X indicate that the 
conclusion concerning the instability of a Skylab-type satellite caused by aero- 
dynamic torque is rather generally valid; neither the specific values of the 
coefficient in the Fourier expansions of the aerodynamic torque nor the specific 
value of the equilibrium angle ^ enter into the final conclusion. Of course, • 
the form of the expansions (i.e., the set of the nonvanishing coefficients) and 
the form of the equilibrium solution do influence the final result. Further- 
more, it is a general experience that perturbations of the considered "pure" 
problem, such as asymmetries of the satellite, a noncircular orbit, the basic 
nonlinearity of the problem, etc., will increase the instability. Thus, we can 
state with reasonable justification that a satellite of Skylab's type is not 
stable under the influence of high altitude aerodynamic torques. 
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